Draft version December 27, 2012 

Preprint typeset using I^-T^X style emulatcapj v. 5/2/11 



(N 

o 

(N 
O 

Q 

w 

Q-i" 

6 

> 



00 

d 

(N 



13 



MAGNETICALLY AND BARYONICALLY DOMINATED PHOTOSPHERIC 
GAMMA-RAY BURST MODEL FITS TO FERMI LAT OBSERVATIONS 

Peter Veres 1 ", Bin-Bin Zhang 2 , Peter Meszaros 1 

(Dated: December 27, 2012) 
Draft version December 27, 2012 

ABSTRACT 

We consider gamma-ray burst models where the radiation is dominated by a photospheric region 
providing the MeV Band spectrum, and an external shock region responsible for the GeV radiation via 
inverse Compton scattering. We parametrize the initial dynamics through an acceleration law r oc , 
with fi between 1/3 and 1 to represent the range between an extreme magnetically dominated and a 
baryonically dominated regime, depending also on the magnetic field configuration. We compare these 
models to several bright Fermi LAT bursts, and show that both the time integrated and the time 
resolved spectra, where available, can be well described by these models. We discuss the parameters 
which result from these fits, and discuss the relative merits and shortcomings of the two models. 



1. INTRODUCTION 

Observations in the GeV range of gamma-ray bursts 
(GRBs) by the Fermi LAT (jAtwood et alJl2009f ) instru- 
ment have uncovered new peculiarities of these objects, 
which were unknown previously. Among the minority of 
bursts which show significant emission above 100 MeV 
(|The Fermi Large Area Telescope Team et al.N2012D . we 
concentrate here on a handful of bursts with substantial 
detected GeV emission to construct models aimed at in- 
terpreting the radiation and the basic parameters of the 
bursts. 

The bursts with observed GeV emission have rather 
diverse spectra. So far, because of relatively low num- 
ber of such GRBs, we restrict ourselves only to bursts 
which have a large enough data coverage to make de- 
tailed model fits meaningful. In general terms, the GeV 
emission in some bursts can be described through an ex- 
trapolation of the lower energy spectrum, without strong 
evidence for departures from this, while in other bursts 
an extra power law component emerges at GeV ener- 
gies, in addition to the spectral component describing 
the MeV radiation. An important new feature present 
in both short and long bursts is the presence of a delay 
between the MeV and the GeV radiation. Furthermore, 
some bursts also exhibit strong spectral evolution. 

The origin of the GeV component has been attributed 
to leptonic mechanisms in the external forward shock, 
e.g. in the r adiative regime in a pa ir enriched interstel- 
lar med ium dGhisellini et~aTl 12009ft or in the adiabatic 
regime ([Kumar fc Barniol Duranl 12009ft. or alternatively 
to hadronic cascades (e.g. lAsano et al.ll2009f) , among oth- 
ers. In the external forward shock, inverse Compton 
scatterin g and Klein- Nishina effects were considered by 
iHe et al.l ()2012l ) , and it has been argued that at least the 
first few seconds of the early high energy emission can- 
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not be attribute d to the forward shock (|He et alj|2011t 
iGao et aT1l2009ft . Thus, it is important to consider both 
the prompt (largely MeV) emission and the GeV emis- 
sion jointly, as they interact with each other. Two recent 
developments which are relevant for the prompt emission 
are the increasing consideration of magnetically domi- 
nated models for the GRB outflow, and of photospheric 
models for the formation of the pro mpt radiation spec- 
trum (e.g. IMeszaros fe Gehrelsl 12011 ). 

Previously, we have shown that some of the bright- 
est LAT detected burst s can be modelled by a m agnetic 
two component model (jVeres fc Meszarosl[2012l ) (|VM12l . 
hereafter). We showed that within this model, without 
going into detailed model fits, it is possible to incorporate 
the seemingly different behaviors of the high-energy spec- 
trum (which is either a single Band function or a Band 
function with an extra high energy component) into a 
single framework. We also addressed the non-detection 
problem of bursts by LAT, namely some of the GRBs 
should have been detected by LAT if the spectrum con- 
tinues smoothly from the MeV regime, or generally there 
is lack of LAT detected bursts compared to initial esti- 
mates, based on the simple extrapolation of the MeV 
emission. We have also presented a generalized formal- 
ism for the dynamics of GRBs valid for a generic acce l- 
eration law T cx r^ with 1/3 < fi < 1 (jVeres et al.ll2012ft . 
where fi can be understood as a parameter accounting for 
the bulk Lorentz factor T(r) averaged across the effective 
cross section of the jet at that radius. 

Here we present detailed numerical fits of the MeV and 
GeV data to specific physical models of GRBs based a 
photospheric origin of the prompt emission in two dy- 
namic regimes, i.e. the extreme magnetically dominated 
regime and the usual baryonically dominated regime, 
with inclusion of the inverse scattering effects (both self 
and external) incurred in the external forward and re- 
verse shocks. This is one of the few instances of fitting 
GRB spectra with physically de t ailed models (for oth er 
examples see iRvde et all (|2010D ; iBurgess et ail (|2011ft ). 
We construct the models based on the generalized dy- 
namics described above, and convolve them with the re- 
sponse function of the Fermi detectors, showing that they 
give a good statistical description of the data. 
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In this paper, we analyze four representative bursts, 
which capture the diversity observed in LAT: GRBs 
080916C, 090510, 090902B and 090926A. We model their 
spectrum with the two component magnetic (p = 1/3) 
model and demonstrate that it is a viable model to inter- 
pret the spectra from approximately 10 keV to 10 GeV. 
We also show that the baryonic model (p, = 1) gives 
similarly good fits to the data, as quantified by specific 
goodness of fit criteria. 

We present our basic model in §2. The Fermi obser- 
vational data for the four GRBs and the fitting method 
are described in §3. The fitting results are presented in 
§4. We summarize the results and compare the bary- 
onic with the magnetic model, with a discussion on the 
physical implications, in §5. 

2. PHOTOSPHERIC MODEL WITH ARBITRARY 
ACCELERATION PARAMETER 

Our basic model consists of a dissipativc photosphere, 
whose location depends on whether the dynamics there 
is kinetic energy dominated or magnetically dominated, 
and an external shock region, where reverse shocks may 
be present as well as a forward blast wave. The bulk of 
the prompt MeV radiation arises from the neighborhood 
of the photosphere in the form of a nonthermal compo- 
nent and depending on the model, a separate, subdom- 
inant thermal component can also be evident. At the 
deceleration radius, shock accelerated electrons upscat- 
ter the prompt photons to produce a component peaking 
in the GeV range. 

In the p resent paper we generalize the model presented 
in IVM12I to cover a wider range of dynamic behavior 
stretching from the extreme kinetic energy dominated to 
the extreme magnetic dominated cases. The acceleration 
of the ejecta is described by the expression (TTJ) for the 
bulk Lorentz factor 



r(r) oc 



r M if r < r sat 
const, if r sat < r < r dcc . 



(1) 



Here the exponent ji is taken as /i = 1 /3 when the dy- 
namics is governed by dissipation of the energy stored 
in magnetic fields in the presence of a sub-dominant in- 
ertia of the baryonic co mponent (e.g. iDrenkhahnl l2002t 
IMeszaros fc Reesl 120 111 ) , while if the baryonic compo- 
nent dominates, we have ji = 1. Intermediate values 
of ji can arise in magnetic outflows with certain symme- 
tries, and the exact value of \i also depends on the mag- 
netic field configuration (e.g . iTchekhovskov et all 120101 : 
iMcKinnev fc U zdcnskv l201lD . r sat and rd ec are the sat- 
uration (coasting) and the deceleration radii, and with 
this paramctrization the radius where the ejecta reaches 

the coasting Lorentz factor is r sat = roT}" . For simplicity 
we consider the transition from the acceleration phase to 
coasting as a broken power law functi on, although this 
trans ition occurs more gradually (e.g. IKobavashi et all 
[19991) . 

In what follows is it is useful to give the dependence 
of the physical parameters in terms of a general ji. The 

4 Depending on model details, the spectrum may start form- 
ing below the phot osphere, and further contributions can arise 
from outside of it f Pc'cr et al. 2006; Lazzati & Bcgclmar] 120101 ; 
IBelob orodov 2010; Giannio SI I201U) ; for our purposes here, we can 
assume that the bulk of the radiation arises at the photosphere. 



main difference between the models is whether the photo- 
sphere occurs below or above the saturation radius. This 
is governed by the relation of rj and a quantity tjt given 
by eq.([2|). Here r\T is a limiting value for the Lorentz 
factor, which separates the two cases of rj > r/T (photo- 
sphere below saturation) and r\ <r\T (photosphere above 
saturation, where tjt is a function of the dynamics, de- 
pending on ji. This limiting Lorentz factor is obtained 
from setting r sat = r ph for a general ji, generaliz ing the 
related value introduced in IMeszaros et al.l (|1993| ) , 
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m/(i+3a0 



(2) 



The photosphere occurs in the acceleration phase if rj > 
ifr, typical for a magnetically dominated (/x = 1/3) case, 

where t\t — 150 I^%Tq\' 6 . The photosphere is in the 
coasting phase for rj < tjt which it is typical for baryonic 
cases (n = 1) for which t\t — 1900 L^r The two 
radii are equal for ji « 0.6 if other nominal parameters 
are left unchanged. The photosphere will occur at 
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VT /fi (VT/v) 1/{1+2fi) if r)>r)T 
Vt^^t/v) 3 if V <Vt 



(3) 
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noting that F p h carries an r ph dependence for rj > 
TjT- The Lorentz factor at the photosphere is T p h = 

(Wro)" « LM/(2M+l) J? -/ 1 /(2^+l) r - M/(2M+1) for ^ > VT 

and r p h ~ fj in other cases. The shell has initially 
a constant width, Ar = ct w sa 3 x 10 t\ cm, where 
t w is the time during which the central engine is ac- 
tive. At a later stage the shell starts to spread radially 
when the when the relative width 8r/r ~ Av/c ~ l/2r 2 
becomes of the same order as the initial shell width 
Ar/r. Depending on the regime, this can happen at 

fsproad = ( 2ct tu /r M J if r < r sat and 2rj 2 ct w \i r > 

r sat . For most of the parameters used here the spread- 
ing occurs well beyond the saturation radius, r sprea( j = 
6 x 10 17 ?7g OO £i cm, and it can be neglected close to the 
photosphere. The scaling of the baryon number density 
in the comoving frame is 



n' b {r) = L/4Trr 2 m p c 3 T(r)rj. 



(5) 



This results in a dependence: n' b (r) oc V (r) oc r 2 M 
if r < ?' sat or r~ 2 if r > r sat . 

Prompt spectrum: Observationally the prompt cmis- 
sion spectrum is described by a Band function 
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of the peak of the spectrum may be due to effects related 
to electron heating by nuclear collisions, reconnection, or 
scmirclativistic shocks, all of which give the peak energy 
at the right order of magnitude. A detailed radiative 
transfer would be model-dependent and too cumbersome 
for performing numerical fits to the data. For the pur- 
poses of this study we choose a simple prescription where 
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the peak energy is determined by the synchrotron peak 
from mildly relativistic shocks around the photosphere^. 

It should be noted however, that the above is just 
a simple heuristic prescription for the peak, the shape 
being taken to be the Band shape. The Band spec- 
tra and peaks deduced in the previously mentioned nu- 
merical simulations arise from distortions of the original 
black-body component formed around the launching ra- 
dius through Compton scattering and synchrotron radia- 
tion, which escapes at the photosphere. In these cases no 
additional black-body component is expected below the 
peak. In our model considered here the peak is ascribed 
to synchrotron radiation from electrons accelerated near 
the photosphere, which introduces a dependence on the 
magnetic field; and the advected blackbody component is 
also released, but at energies typically below the peak, for 
magnetically dominated dynamics. There are possible 
shortfalls for a pure synchrotron interpretation, such as 
the low energy spectral slope observations which show in- 
consistencies with the synchrotron model when electron s 
are in the fast cooling regime (e.g. iPreece et al.lll998[ ). 
While recent analyses (jSakamoto et al.l 120 1 If ) show that 
the fraction of bursts showing problems is smaller then 
previously thought, there is still a non- negligible frac- 
tion of bursts where other mechanisms may need to be 
invoked. This could be the Compton upscattcring men- 
tioned above, or modifications of the simple synchrotro n 
mechanism, e . g. su ch as discussed bv iMedvedevi (|2006f ): 
iDaigne et al.l (j 2 1 If ) . A first principles discussion of such 
cases for our fitting exercise undertaken here on individ- 
ual bursts would be too extensive, and for this reason we 
have adopted the simple prescription mentioned above. 

The peak of the prompt spectrum is then given by £f, 

(3 ge B ph /4™ eC ) 7o 2 ph [r ph /(l + z)], where B' " 
(taking n' b from Eq. [5]) and J e ,ph oc IV 



if ''sat < r < r S p rea d. This results in a dependence of the 
observed temperature as: 



oc n b I r 



£b oc 



Ivwtj ^+2 r ^+ 2 T3/(l + z)ii n > r\ T 



L- 1 /vr3/(i + z) 



if Tj < TjT- 



(6) 



Here we have assumed an efficiency of 0.5 (jGianniosI 
l201ll) . The two cases imply widely different values for 
the same set of parameters. This is mainly caused by 
the strong dependence on rj when the photosphere arises 
in the coasting phase (rj < t]t)- Note furthermore, that 
the photosphere loses any memory of the acceleration in 
the same case and no [i dependence is present. Also in 
this case, the peak energy is independent of the launching 
radius. We take the low and high energy photon number 
indices from observations (a ~ 1, /3 « 2.2), and in the 
fitting sections we directly fit for these. 

Thermal component: Even though the main emis- 
sion component is nonthermal, one expects a thermal 
component to arise from the photosphere. There are 
hints of this separa te thermal component emerging be- 
low the Band peak (IPage et al.ll2011l:lGuiriec et a.l.ll2011h 
iMcGlvnn fc Fermi Collaboration! 120121 ). The comoving 

temperature goes as T" oc r if r < r sa t and r < 3 

5 Shocks are expected as result of magnetic reconncction, as 
well as possibly from outflow Lorentz factor fluctuations, the par- 
ticle acceleration time being of the sam e orde r for reconncction 
and Fermi acceleration Mcszaros & Recs (20fl), and the spectral 
peaks are in both cases in the MeV range, as for baryonic thermal 
photosphere models. 
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,L- 5 / 1 V /3 ro /6 /(l + . 
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External shock-related components: The ejecta starts 
to decelerate when it has plowed up approximately 1/r/ 
times its own mass from the interstellar matter. This oc- 
curs at the deceleration radius rd cc = 4.8 x l0 16 L t ' 53 (1 — 

C r ) 1 / 3 t} / 3 3 n " 1 ^ 3 ?7g Q' /3 cm. Here £ r is the fraction of the 
luminosity in radiation, which for simplicity is fixed to 
0.5 throughout the rest of this article. At the decel- 
eration radius a forward shock develops (and possibly 
a reverse shock as well). The prompt emission is up- 
scattered on the shock-accelerated electrons, giving rise, 
to a > 1 GeV peaked radiation. This component dom- 
inates at sufficiently high energies, but radiation from 
synchrotron self Compton (SSC) from the forward shock 
(FS) may also be relevant here. The deceleration radius 
is larger than the photosphcric and the saturation radius 
for every realistic set of parameters. 

3. DATA AND FITTING METHOD 
3.1. Data 



Currently there are at least 
energies 



40 LAT- 

detected GRBs at energies > 100 MeV 
(lOmodei fc the Fermi LAT collaboration! 120111: 
IZheng et al.1 I2012T) . Out of these, there are only a 
handful of bursts with enough data at GeV energies to 
derive meaningful constraints on the model parameters 
(jZhang et al.l 12011! ) . In this paper we choose four of 
such bright LAT GRBs, namely, GRB 080916C, 090510, 
090902B and 090926A as our sample. 

We processed the Fermi GBM and LAT data usin g 
the same method as presented in iZhang et al.l (|2011l) . 
For each burst, we choose the 1-2 brightest GBM spec- 
tra and 1 or 2 BGO spectra together with LAT spec- 
trum to compare with our model. We have fixed the 
BGO normalization factor, accounting for differences be- 
tween detectors' effective area, to the u sual value of 0.9 
(|Abdo fc the Fermi collaborationll2009l) . 

By varying the free parameters, we get the best models 
which can reproduce the observed count spectrum. For 
the total duration of the bursts we managed to fit GRBs 
080916C and 090510 with our model. For GRBs 090926A 
and 090902B, our model did not give a satisfactory fit 
when considering the entire duration of the burst, but 
time resolved spectra were well fitted by our model. We 
put this on the account of the strong spectral evolution 
present in these two GRBs as published in their initial 
discovery papers. 

3.2. Fitting method 

We have taken the theoretical model described above 
with 9 free parameters. We fit the data using the sim- 
plex method imp lemented in the amoeba routine (see 
IPress et al.l I2007D modified, in some of the cases with 
simulated annealing (Erik Rosolowsky, private commu- 
nication). This addition makes the method less sensitive 
to local minima. The parameters describe the spectrum 
of a photosphcric model (cither magnetic or baryonic). 
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TABLE 1 

Model parameters for the time-integrated fits with a 
magnetically and a baryonically dominated model for two 
bursts. Stars mark parameters for which we could not 
determine the errors. 



The procedure is to convolve the model with the de- 
tectors' response functions and compare this result with 
the observed counts. For this purpose, we have written 
our own IDL codes to handle both the observed count 
spectrum, the instrumental response matrix (DRM) and 
read in the multi-parameter theoretical model to carry 
out the fitting. 

Because the number of fitting parameters is large, 
we introduce some additional constraints, in accordance 
with the model. We require the deceleration time to be 
of order of few seconds, which corresponds to the delay 
between the MeV and GeV components. This is achieved 
by introducing a penalty in the goodness of fit for t& cc > 
few s. We vary one by one each of the parameters around 
the best fit to calculate their errors. When the \ 2 statis- 
tic increases by the amount characteristic for the one- 
sigma limit for given degrees of freedom, we fix the range 
and thus the errors for that parameter. 
Parameters constrained: In a situation when the spec- 
trum is described by either an extremal magnetic (p = 
1/3) or baryonic (/i = 1) model, the varying parame- 
ters we set out to constrain are: the break energy (eb, 
see Eq. of the prompt emission, the low- and high 
energy indices (a and (3) of a Band function (N e (e) oc 

e -a e -e(2- a )/e b for £ < £& an diV £ ( £ ) (X if £ > S b ), the 

relative Lorentz factor of the reconnecting regions in the 
photosphere F r (> 1), the total luminosity L t , interstellar 
density (n e ), index of the electron power-law distribution 
(p, where N e {l) oc 7~ p ) and the fraction of the energy 
in the magnetic field and electrons at the deceleration 
radius (cb,fs and e e ). We do not to include variations 
in e.g. magnetic field in FS and RS because it would 
not be possible to constrain them and we assume that 
£b,fs = £ b.rs- The components included in this work 
are: prompt Band component (PR), blackbody compo- 
nent (BB), forward- and reverse shock synchrotron (FS, 
RS), forward and reverse shock synchrotron self Comp- 
ton (FS-SSC, RS-SSC), external inverse Compton of the 
prompt photons on the forward and the reverse shock 
(FSEIC, RSEIC). 

The prompt emission is represented Band func- 
tion, the forward and reverse shock synchrotron and in- 
verse Compton components as multiple broken power 
laws. We have applied a phcnomenological smoothing at 
the breaks of the type ((x / x hrcak ) llb + {x/x hYcak )^ b ) 1 / b . 



Here b is an arbitrary constant describing the smooth- 
ness of the break (we set b = —4), ii j2 are the power law 
indices and the break occurs at screak- The prompt emis- 
sion may also have a high energy spectral cut-off, at an 
energy depending on the acceleration mechanism at the 
photosphere, and whether or not pair creation is present. 
Here we have taken a phenomenological cut- off value at 
appro ximately r p hTO e c 2 /(l + z) (as also in iThompsonl 
(|1994l) ). Both without pair formation and for pair for- 
mation estimated from the photons above threshold es- 
timated from an extrapolation of the observed spectral 
slopes in the bursts considered here, a cutoff close to the 
above value is obtained. The energy coverage is generally 
sparse at a putative prompt cutoff energy in the range 
of tens of MeV, making it hard to discern whether there 
is a pair cutoff or a cutoff due to the particle accelera- 
tion mechanism. For simplicity, here we use the model 
without pair creation for the purposes of the numerical 
calculation. This cutoff can constrain the photospheric 
Lorentz factor (r p h). Concerning the GeV component, 
we identify the high energy cutoff (and peak) with the 
Klcin-Nishina cutoff. 

The strength of the GeV component is significantly 
influenced by the amplitude of the prompt component 
(responsible for the intensity of seed photons for the EIC 
scattering), e e and n ext , while cb,fs plays a role in some 
of the cases. The shape of the GeV emission is defined 
by a combination of a and p, which makes it possible to 
constrain these parameters. The peak energy of the GeV 
emission depends on 77, e e and p. The coasting Lorentz 
factor, r\ can be constrained by the prompt emission cut- 
off at tens of Me Vs. 

The model in general has difficulty in constraining 
£b,fS: as this parameter affects only the low-energy 
(~ 10 kcV) behavior of the spectrum through the FS 
synchrotron (and to a lesser extent by RS-SSC emis- 
sion). As a consequence this energy range can be used 
to constrain the parameter, but in many cases, however, 
there are no signs of an extra component emerging at 
low energies. We note that, considering the small num- 
ber of high-energy photons, the fact that many of the 
combined GBM-LAT spectra are consistent with a sin- 
gle Band function does not necessarily exclude additional 
components. The two components (photospheric and ex- 
ternal shock-related) can align to form a continuous com- 
ponent without excessive fine tuning. This requirement 
is relaxed when we need to fit bursts with additional 
components. 

4. INDIVIDUAL GRBS 

4.1. GRB 090510 fit for time integrated spectrum 

This is a short burst (jAbdo et al.ll2009D with GeV emis- 
sion. According to our modelling the time integrated 
spectrum can be described by a Band function and the 
high energy photons are dominated by the RSEIC emis- 
sion. The low-energy excess is best modelled by a FS 
component in the magnetic case and by a blackbody com- 
ponent in the baryonic case. See Figs. [I] and [3] 

4.2. GRB 080916C fit for time integrated spectrum 

GRB 08 0916C is one of the most intensive burs ts ever 
detected ()Abdo k the Fermi Collaboration! 120091) . Its 
spectrum is consistent with the Band function through- 
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Fig. 1.— GRB 090510 fit with a magnetic model 
(x 2 /dof=107.6/97). In the upper figure lines represent the 
model. The points with the error bars can be thought of as 
measured values, though they are themselves model dependent 
(i.e. for different models they are not the same). They are 
calculated from the model convolved with the detectors' responses 
(step curve on lower figure) and their relation to the observed 
counts (crosses on lower figure). For brevity we do not show the 
plot with the convolved spectra for other bursts. 



out the Fermi energy range and in all the analyzed inter- 
vals ( though there was weak evidence of an ex tra compo- 
nent (jAbdo fc the Fermi Collaboration! 120091 ) ). The ex- 
tremal magnetic and baryonic fits are shown in Fig. |3l 



GRB 080916C 



E 



10"" 




t 








L=7.6E+54 erg/s _ 

e =5.6E-02 
■e r =3.1E-02 

f=C50 : 
p=2.B7 - 
o=0.5 














r,=i .3 _ 


10~ 7 


: \>i 

\ ' / 

V. 

/ \ 








\ 
\\ 
\\ 
\ v 

/ 1 
/ 1 






10~ s 


1 

1. ,\ . / 


' \ 







10° 10 2 10" 10 6 10 a 

e [keV] 
GRB 08091 BC 




10° 10 J 10" 10 6 10' 

, [k=V] 



Fig. 3.— Fit for GRB 080916C, for the entire interval, with a 
magnetic model (above, x 2 = 113.6, d.o.f.=98) and baryonic model 
(below, x 2 = 129.7, d.o.f.=98). The figures show the models and 
the data points obtained by scaling the measured count rates with 
respect to the convolved models. 
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Fig. 2.— GRB 090510 fit with a baryonic model 
X 2 /dof= 107. 1/97). Other notations are the same as Fig. 



4.3. GRB 090926 'A fit for the time resolved spectrum 

lAckermann fc the Fermi co llaborat ion! ()2011[ ) fitted 
the spectrum of this burst with a Band function and an 
extr a power-law with a high energy cutoff. Using data 
from iZhang et al.l (|2011f ) the whole interval of the burst 
cannot be well fit with a simple analytical model, sug- 
gesting strong spectral evolution. We thus use the time 
resolved data for this GRB. We divided the burst in 4 
intervals (with limits at 0, 4, 7, 11.5 and 19 s with re- 
spect to the trigger time) . The first interval has no GeV 
range emission, and we drop it from our analysis. 

We found that both the magnetic and the baryonic 
photosphere model can give an accurate description of 
the resolved data (see Fig. 0] for a superposition of mag- 
netic and baryonic models and the models and data for 
the first interval on Fig. [5]). 

Reduced \r values are in the range 1.4 — 1.8 for the 
magnetic model, and 1.2 — 2.5 in the baryonic case. Gen- 
erally, e e , €b and terminal Lorentz factors (rj) are system- 
atically lower in the baryonic case, the external density 
is systematically higher, while the other parameters are 
roughly in the same range for both cases. We note as a 
general trend, the magnetic model yielded better fits for 
the data. 



4.4. GRB 090902B fit for the time resolved spectra 
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Fig. 4.— GRB 090926 A with three intervals fitted using the 
magnetic (top) and baryonic (bottom) photosphere model. We 
have displayed only a small number of components. The numbers 
indicate the intervals: 1: 4-7 s, 2: 7-11 and 3: 11.5-19 s. The 0-4 s 
period could not be fitted, because no LAT photons were detected. 
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Fig. 5. — Examples of magnetic (above, x 2 = 182.9, d.o.f.=115) 
and baryonic (below, x 2 = 184.5, d.o.f.=115) models with the data 
overplotted for GRB 090926A's 4-7 s interval (see also Fig. gj). 



GRB 090902B (jAbdo fc the Fermi collaboratioij[2009l) 
has a relatively sharp Band spectrum with an underlying 
power law component. While a simple Planck function is 
too narrow to fit the spectrum, a multicolor blackbody 
and an underlying po wer law can fit the data equally 
well (jRvde et al.H2010t) . The fit for the entire duration of 
the burst did not give an acceptable fit for neither of the 
models, possibly due to the strong spectral evolution. It 
has been show n that this burst in volves a strong spectral 
evolution (e.g. lZhang et alJl2011l ) 

In order to account for spectral evolution, 
we have defined 10 intervals bracketed by the 
0,5,7.3,9.7,12,14.3,16.7,19,21.3,23.7 and 26 s marks 
relative to trigger (see Fig. [6]). 
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Fig. 6. — GRB 090902B magnetic models (upper) and baryonic 
(lower). The time intervals depicted: 0-5 s: gray-dotted, 7.3-9.7 
s: black-long dashed, 9.7-12 s: gray-dash-dottcd, 12-14.3 s: gray- 
dash-dot-dot-dot, 14.3-16.7 s: black-continuous, 16.7-19 s: gray- 
continuous, 19-21.3 s: black-dashed, 21.3-23.7 s: light gray-dash- 
dotted. 

e e varies in the range « 0.001 — 0.7, and it is similar 
for both models. The terminal Lorentz factor is system- 
atically higher in the magnetic case and it is above 400 
in all intervals, while in the baryonic case, it is below 
250. The external density, is systematically lower in the 
magnetic case, while the magnetic parameter at the de- 
celeration radius and the semirelativistic shocks' Lorentz 
factor in the photosphere are systematically higher in the 
magnetic case. These parameters are all within the nor- 
mal ranges found previously from fitting afterglow and 
prompt emission data on a wider sample of bursts. Here, 
we do not find any obvious trend for an evolution of the 
parameters throughout this burst. For an example fit for 
this burst's first interval, see Fig. [7] 

5. DISCUSSION 

Possible relation between fi and observed trends.- In 
our fits we have purposefully avoided introducing any of 
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Fig. 7. — Examples of magnetic (above, \ 2 = 102.3, d.o.f.=92) 
and baryonic (below, \ 2 = 96.7, d.o.f.=92) models with the data 
overplotted for GRB 090909B's 0-5 s interval (see also Fig. \E\- 
gray-dotted curve). 



the observed phenomenological correlations claimed for 
GRBs, since these have errors and biases which require 
extensive consideration and would complicate the analy- 
sis. However, it is interesting to speculate on the possible 
values of fi in the context of some of these phenomeno- 
logical correlations. 

In the J] > tjt case it seems simple considerations do 
not yield a consistent value for fi, and this needs to be 
examined further. In the other case (r) < ?y<r), there is 
no fi dependence, but we only need to assume one of the 
relations to obtain a pproximately th e other one; e.g. if 
we assume r\ cx L 3 (ILii et all 120121) . we get Eb oc L 0A , 
fairly close to the Yonetoku relation. 

Although we did not find any cases where the main 
MeV episode was provided by the thermal component, 
such cases have been discussed based based on more elab- 
orate r adiation tran s fer m o dels, e.g. by iPe'er fc Rvdd 
(poll : iBeioborodovl (pOll : IVurm et all poll " and oth- 
ers, typically for baryonic dynamics. In our simplified 
treatment the thermal component was taken as a sim- 
ple blackbody function, instead of a broa der multicolor 
blackbody as in, e.g. (|Rvde et al.llMof ). Such multi- 
color fits including our other variables would be much 
more cumbersome, but in principle, if we make a quick 
calculation similar to that in the previous paragraph, us- 
ing Eq. [7] instead of Eq. |H1 this would imply /i ps 0.5. 

It is interesting to note that for the magnetically dom- 
inated model {rj > t\t and /i = 1/3), the peak energy 
does not depend on either the luminosity and the termi- 
nal Lorentz factor, and has a value of ~ 300 keV, cor- 
responding to the c lustering, observed by e.g. BATSE 
(|Kaneko et al. 1 120061) . There is a weak dependence on r 
and redshift. 

Comparison of the baryonic photosphere model with the 



magnetic model.- The main difference between the mag- 
netic and the baryonic model in terms of the physical 
parameters involved in the GRBs is the initial depen- 
dence of the Lorentz factor on the radius (see equ. [TJ. 

The main manifestation of this difference in terms of 
the observed spectrum is the strength of the thermal 
component from the photosphere. Because of the dif- 
ferent dynamics, the magnetic model results in a sub- 
dominant thermal component, while in the case of the 
baryonic dynamics the black body component can have 
a more prominent contribution. 

The main goal of this article has been to present a gen- 
eralized formalism for models with an arbitrary acceler- 
ation of the GRB ejecta, valid for magnetic and bary- 
onic dynamics. By means of an analysis of the brightest 
bursts, and properly accounting for the instrumental ef- 
fects, we have shown that both an extremal magnetic 
and a baryonic model can provide adequate fits to the 
spectrum. Although in some of the cases the magnetic 
model provides a marginally better fit to the data, the 
difference is not significant enough to decide in favor of 
either extremal model. 

Constraints on the parameters.- As we can see from 
Table 1. we can constrain some of the physical parame- 
ters of the models, except for the magnetic parameter at 
the deceleration radius, and in some cases the interstellar 
density is also difficult to constrain. 

We found that the burst spectra are consistent with 
both the magnetic and the baryonic models presented 
here. A possible discrimination would require further 
theoretical considerations, or else other observational 
considerations; e.g. the length of the deceleration time 
should be the same order as the observed delay of the 
GeV emission. Another discriminant would be how well 
a model can reproduce the observed correlations between 
the rest-frame peak energy and e.g. luminosity would. 
As we indicated above, the extreme magnetic photo- 
sphere model appears to be better in accommodating 
these constraints. 

One of the main differences between the two models is 
the systematically higher value of the coasting Lorentz 
factor (LF) in the magnetic model. This can be under- 
stood as follows: by applying two models to the same 
data the photospheric LF is the best constrained by the 
observations. While in the baryonic case this is just the 
terminal LF, in the magnetic case the ejecta is still ac- 
celerating at the photosphere, having a larger coasting 
LF. 

In GRB 090510 lAbdo et al.l (|2009[ ) noted a correlation 
between the keV range radiation and the GeV radia- 
tion (but no correlation with the MeV radiation). Based 
on their spectrum, the additional power-law would be 
the responsible for explaining this correlation. In the 
framework of this generalized photospheric model, such 
a correlation can be explained if these two components 
have a common source: e.g. the keV emission is as- 
cribed to the forward shock (FS) synchrotron and the 
GeV to the forward shock synchrotron-self Compton (FS- 
SSC); or if a reverse shock is present (as proposed in 
this work), the keV radiation can be due to the reverse 
shock synchrotron-self Compton (RS-SSC) and the GeV 
is produced by the reverse shock external inverse Comp- 
ton (RS-EIC). 

Conclusions.- In summary, we have shown that a sim- 
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pic dissipativc photosphere model of the prompt emission 
of GRB based on synchrotron radiation can fit well the 
joint Fermi GBM and LAT data for the four brightest 
bursts with sufficiently detailed information. We have 
introduced a general formalism for treating the jet accel- 
eration dynamics ranging between an extreme magneti- 
cally dominated to a baryon dominated dynamics regime. 
Both types of dynamics give acceptable fits to the data 
on these four bursts. Additional considerations involv- 



ing statistical trends over a much larger sample of bursts 
might favor slightly the magnetic models, but firm con- 
clusions would require detailed fits over much larger sam- 
ples than undertaken here. 

We acknowledge partial support from NASA 
NNX09AL40G and OTKA K077795 (PV, PM) and 
NASA SAO SV4-74018 and NASA SAO GO1-12102X 
(BBZ). We thank Kazumi Kashiyama, Yi-Zhong Fan and 
the anonymous referee for comments and suggestions. 



REFERENCES 



Abdo, A., & the Fermi Collaboration. 2009, Science, 0036, 8075 
Abdo, A. A. et al. 2009, Nature, 462, 331 

Abdo, A. A., & the Fermi collaboration. 2009, ApJ, 706, L138, 
0909.2470 

Ackermann, M., & the Fermi collaboration. 2011, ApJ, 729, 114, 
1101.2082 

Asano, K., Guiriec, S., & Meszaros, P. 2009, ApJ, 705, L191, 
0909.0306 

Atwood, W. B. ct al. 2009, ApJ, 697, 1071, 0902.1089 
Band, D. et al. 1993, ApJ, 413, 281 

Bcloborodov, A. M. 2010, MNRAS, 407, 1033, 0907.0732 

. 2012, ArXiv c-prints, 1207.2707 

Burgess, J. M. et al. 2011, ApJ, 741, 24, 1107.6024 
Daigne, F., Bosnjak, Z., & Dubus, G. 2011, A&A, 526, A110, 
1009.2636 

Drenkhahn, G. 2002, A&A, 387, 714, |arXiv:astro-ph/0112509 



Gao, W.-H., Mao, J., Xu, D., & Fan, Y.-Z. 2009, ApJ, 706, L33, 
0908.3975 

Ghiscllini, G., Ghirlanda, G., Nava, L., & Celotti, A. 2009, ArXiv 

c-prints, 0910.2459 
Giannios, D. 2011, ArXiv c-prints, 1111.4258 
Guiriec, S. et al. 2011, ApJ, 727, L33, 1010.4601 
He, H.-N., Liu, R.-Y., Wang, X.-Y., Nagataki, S., Murase, K., & 

Dai, Z.-G. 2012, ApJ, 752, 29, 1204.0857 
He, H.-N., Wu, X.-F., Toma, K., Wang, X.-Y., & Meszaros, P. 

2011, ApJ, 733, 22, 1009.1432 
Kaneko, Y., Preece, R. D., Briggs, M. S., Paciesas, W. S., 

Meegan, C. A., fe Band, D. L. 2006, ApJS, 166, 298, 

|arXiv:astro-ph/0601188| 
Kobayash i, S., Piran, T~~fc Sari, R. 1999, ApJ, 513, 669, 

|arXiv:astro-ph/9803217| 
Kumar, P., fc Barniol Duran, R. 2009, MNRAS, 400, L75, 

0905.2417 

Lazzati, D., & Begelman, M. C. 2010, ApJ, 725, 1137, 1005.4704 
Lii, J., Zou, Y.-C, Lei, W.-H., Zhang, B., Wu, Q., Wang, D.-X., 

Liang, E.-W., & Lii, H.-J. 2012, ApJ, 751, 49, 1109.3757 
McGlynn, S., & Fermi Collaboration. 2012, in Gamma-Ray 

Bursts 2012 Conference, ed. A. Rau and J. Greiner, 012 
McKinney, J. C, & Uzdensky, D. A. 2011, MNRAS, 1766, 

1011.1904 



Medvedev, M. V. 2006, ApJ, 637, 869, |arXiv:astro-ph/051 0472 
Meszaros, P., & Gehrels, N. 2012, Research in Astronomy and 

Astrophysics, 12, 1139, 1209.1132 
Meszaros, P., Laguna, P., fe Rees, M. J. 1993, ApJ, 415, 181, 

|arXiv:astro-ph/9301007| 
Meszaros, P., & Rees, M. J. 2011, ApJ, 733, L40+, 1104.5025 
Omodei, N., & the Fermi LAT collaboration. 2011, Talk at Fermi 

Meeting, Stanford U. 
Page, K. L. et al. 2011, MNRAS, 416, 2078 
P e'er, A., Meszaros, P., fc Rees, M. J. 2006, ApJ, 642, 995, 

|arXiv:astro-ph/0510114| 
Pe'er, A., fc Ryde, F. !201l, ApJ, 732, 49 
Pe'er, A., Zhang, B.-B., Ryde, F., McGlynn, S., Zhang, B., 

Preece, R. D., & Kouveliotou, C. 2012, MNRAS, 420, 468, 

1007.2228 

Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., 
Paciesas, W. S., fc Band . D. L. 1998, ApJ, 506, L23, 
|arXiv:astro-ph/9808184| 
Press, W. H., Teukolsky, S. A., Vcttcrling, W. T., & Flannery, 
B. P. 2007, Numerical Recipes 3rd Edition: The Art of 
Scientific Computing, 3rd edn. (New York, NY, USA: 
Cambridge University Press) 
Ryde, F. ct al. 2010, ApJ, 709, L172, 0911.2025 
Sakamoto, T. et al. 2011, ApJS, 195, 2, 1104.4689 
Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2010, 15, 749, 
0909.0011 

The Fermi Large Area Telescope Team et al. 2012, ApJ, 754, 121 
Thompson, C. 1994, MNRAS, 270, 480 
Veres, P., & Meszaros, P. 2012, ApJ, 755, 12, 1202.2821 
Veres, P., Zhang, B.-B., & Meszaros, P. 2012, ArXiv e-prints, 
1208.1790 

Vurm, I., Beloborodov, A. M., & Poutanen, J. 2011, ApJ, 738, 77, 
1104.0394 

Vurm, I., Lyubarsky, Y., & Piran, T. 2012, ArXiv e-prints, 
1209.0763 

Zhang, B.-B. et al. 2011, ApJ, 730, 141, 1009.3338 
Zheng, W., Akerlof, C. W., Pandcy, S. B., McKay, T. A., Zhang, 
B., Zhang, B., & Sakamoto, T. 2012, ArXiv c-prints, 1203.5113 



